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ABSTRACT 

We investigate the difference between the host galaxy properties of core-collapse supernovae 
and long-duration gamma-ray bursts (LGRBs), and quantify a possible metallicity dependence 
of the efficiency of producing LGRBs. We use a sample of 16 CC SNe and 16 LGRBs from 
iFruchter et alJ (120061) which have similar redshift distributions to eliminate galaxy evolution 
biases. We make a forward prediction of their host galaxy luminosity distributions from the 
overall cosmic metallicity distribution of star formation. The latter is based on luminosity 
functions, star-formation rates (SFR) and luminosity-metallicity (L-Z) relations of galaxies. 
This approach is supported by the finding that LGRB hosts follow the L-Z-relations of star- 
forming galaxies. We then compare predictions for metallicity-dependent event efficiencies 
with the observed host data. We find that UV-based SFR estimates predict the hosts distri- 
bution of CC SNe perfectly well in metallicity-independent form. In contrast, LGRB hosts 
are on average fainter by one magnitude, almost as faint as the Large Magellanic Cloud. As- 
suming this is a metallicity effect, the present data are insufficient to discriminate between a 
sharp cutoff and a soft decrease in efficiency towards higher metallicity. For a sharp cut-off, 
however, we find a best value for the cutoff metallicity, as reflected in the oxygen abundance, 
12 + log(Q/H)i, ~ 8.7 ± 0.3 at 95% confi dence including systematic uncertainties on 
the calibration of Kobulnick v & Kewlevl d2004l) . This value is somewhat lower than the tra- 
ditionally quoted value for the Sun, but is comparable to the revised solar oxygen abundance 
(Asplund. Grevesse & Sauvall2005l) . LGRB models that require sharp metallicity cutoffs well 
below ~ 1/2 the revised solar metallicity appear to be effectively ruled out, as they would 
require fainter LGRB hosts than are observed. We also discuss the likely implications of the 
still ongoing metallicity 'calibration debate'. 
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1 INTRODUCTION 

Since the first optical afterglow to a gamma-ray burst was 
identified dvan Paradiis et alj|l99l GRB 970228), the restframe 
properties of about 50 host galaxies of long-duration gamma- 
ray bursts (LGRBs) have been identified. Already with a small 
sample it became cle ar that the hosts are preferentially faint 
blue irregular galaxies iFruchter et all 1999; Le Floc'h etld]2003t 
IChristensen. Hiorth & Gorosabel 2004) . This trend was recently 
confirmed by a comparison between hosts of LGRBs and hosts of 
core-collapse supernovae (CC SNe). CC SNe were found in late- 
type galaxies of all morphologies and luminosities, in contrast to 
LGRB hosts, which were st ill found to be almos t exclusively irreg- 
ular and on average fainter dFruchter et alf2 006, hereafter F06). 

Since both CC SNe and LGRBs are considered explosions 
of massive young stars with ages below ~ 30 Myr and proba- 
bly below the typical timescale of star formation in most galax- 
ies, there is apriori no reason to expect a considerable difference 
in their host galaxies. CC SNe should be direct and unbiased trac- 



ers of very recent star formation, although complications may arise 
from dust extinction intrinsic to the host. Furthermore, we expect 
a dependence of the mass threshold for core collapse/black-hole 
formation on metallicity faan. P odsiadlowski & Eg gletonl 1 19951 : 
iHeger et al. 2003; Eldridge & Tout 2004) a nd binary evolution ef- 
fects llBrown. Lee & Bethel 1 19991 : [Brown et alJ l200ll : IPfahl et al.1 
l2002l:|Podsiadlowski et alj|2004bh . 

The suggested explanation for the host preference of LGRBs 
is that they arise from low-metallicity stellar populations which are 
predominantly found in low-luminosi ty galaxies according to the 
luminosity-metallicity relation fe.g [Kobulnickv ^ Kewlevl 1 20041 
hereafter KK04). Also, in se veral progenitor models of GR Bs (in 
particular 'collapsar' models; MacFavden & Wooslev 1999), wind 
mass loss plays an important role in slowing down the rotation of 
massive stellar cores. Since lower metallicity leads to weaker stel- 
lar winds and hence less angular-momentum loss, the stars are more 
likely to retain rapidly rotating cores at the time of the explosion, 
as required in the collapsar model. 
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LGRBs discovered by the Swift mission have a median red- 
shift of 2.75 and apparently hav e a distribution expect ed for an un- 
biased tracer of star formation I Jako bsson et al.ll2006l) . If both the 
cosmic star formation history and the cosmic evolution of metal- 
licity in cold star-forming gas were known with high accuracy, this 
redshift distribution could provide clues about the metallicity de- 
pendence of the LGRB rate. However, the LGRBs discovered by 
BeppoSAX have a median redshift near z ~ 1, and predictions 
for the completeness function are very difficult to make, owing to 
the complexity of GRB trigger algorithms and detector intricacies. 
Even with a much larger sample, the redshift distribution of GRBs 
might be affected by systematic uncertainties such that no strong 
constraints on metallicity can be found. 

Instead, first spectroscopic measurements of the metallicity 
of LGRB hosts at redshifts 0.4 < z < 1 fr om the GHost Stud- 
ies dSavaglio. Glazebrook & LeB orgne 2006) indicate that LGRB 
hosts follow the luminosity-metallicity relation of star-forming 
galaxies. A preference for low metallicity would then require 
LGRBs to favour less luminous galaxies more than general star for- 
mation does. When working with gas-phase metallicities of galax- 
ies, we need to be aware of the calibration debate between two ba- 
sic methods which result in numbers that differ by a factor of ~ 2: 
spectra of low signal-to-noise only facilitate an indirect strong-line 
method (Pagel's R23 indicator, 1979) to measure metallicity, which 
is calibrated against photoionization models. In contrast, high-S/N 
spectra allow for a direct measurement of metallicity taking into 
account t he electron temperatur e measured from weak auroral lines 
(see e.g. iKennicu tt et al. 2003). The T e -based methods find typi- 
cally 0.2-0.5 dex lower metallicities, but have been critized to be 
biased by not taking inhomogeneities in the temperature into ac- 
count. 

In this paper, we wish to derive constraints on the metallicity 
dependence of LGRB progenitors by comparing forward predic- 
tions of host luminosity distributions with the observed host sam- 
ple. We apply a minor modification to the F06 samples by restrict- 
ing them to identical redshift ranges 2 = [0.2, 1] (see Sect. 2). 
Our quantitative predictions of the host luminosity distribution is 
based on luminosity functions and star-formation rates of galaxies 
and their luminosity-metallicity relations (see Sect. 3). Here, we 
need to make sure that all ingredients are on the exact same metal- 
licity calibrator and opt for the strong-line method on the KK04 
calibrator for consistency. We then fold in different prescriptions 
for the metallicity dependence of the LGRB event rate and con- 
strain these by comparing our predictions to the observed sample 
in Sect. 4. In Sect. 5 we discuss the results and the implications of 
the calibration debate. Throughout the paper, we use Vega magni- 
tudes and the cosmological parameters (fi m , Q\) — (0.3, 0.7) and 
Ho = h.70 x 70 km/(s Mpc). We use t he revised solar oxygen abun- 
dance of 12 + log (Q/H) = 8.66 jAsplund. Grevesse & Sauvail 
l2005l see also All ende Prieto. Lambert & Aspl und (2001)). 



2 DATA 

This work is based on the host galaxy data of CC SNe and LGRBs 
collected by F06. The CC SNe wer e all discovered by the Hubble 
Higher Redshift Super nova Searc h (Strolger et al. 2004) in the two 
HST GOODS fields dGiavalisco et alj|2004l) . This search was only 
sensitive to CC SNe at z ^ 1 and may already be inc omplete at the 
upper end of this redshift range (Dahlen et al. 2004). The LGRBs 
were taken by F06 from public HST data of host galaxy observa- 
tions, where the list of objects extends to redshifts above 4. 



At redshifts z < 0.2, the HST GOODS survey has too lit- 
tle volume to detect CC SNe efficiently. The same restriction ap- 
plies to samples of LGRBs with typical gamma-ray luminosities. 
However, the GRB list contains a few very nearby objects with 
extremely low gamma-ray luminosities, which may be viewed at 
larger angles or be intrinsically faint events, and possibly an alto- 
gether different k ind of explosion (see the ensuing discussion on 
GRB 060218, e.g. lSoderberg et alj20o"6T) . In order to exclude these 
peculiar GRB events, which we could only see at redshifts z < 0.1, 
we restrict the existing CC SN and LGRB data to a subsample in 
the redshift range 2 = [0.2, 1.0]. We note, that this is similar to but 
more conservative than the comparison considered by F06, which 
includes LGRBs up to 2 = 1.2. 

As a result, we keep the full sample of 16 CC SNe with a 
mean redshift of 0.63, while the LGRB sample is restricted to 16 
events with a mean redshift of 0.69, so they are on average from an 
extremely similar cosmic epoch. 



3 MODELLING THE HOST LUMINOSITY 
DISTRIBUTION 

We wish to model restframe V^-band luminosity distributions for 
the host galaxies of both types of explosive events in order to com- 
pare them with the observed distribution from the restricted F06 
sample. For this purpose we ultimately need a distribution of star- 
formation density over l/-band host luminosity and potentially a 
metallicity-dependent efficiency function. 

3.1 Galaxy luminosity function and star-formation density 

IWolf et alj d2005l hereafter W05) measured the restframe UV 
(280 nm) luminosity function (LF) of galaxies from a sample of 
almost 1500 galaxies at a mean redshift of 2 ~ 0.7 0. Based on mor- 
phological classifica tions from the GEMS survey dRix et al1l2004l : 
ICaldwell et afl l2006), the LFs were also split by galaxy type. The 
mean redshift of these LFs coincides with that of the host galaxy 
samples for the CC SNe and LGRBs considered here. Hence, they 
can serve as an unbiased description of the general galaxy popula- 
tion at the redshift of the events. The LF <^>(I/28o), is parameterized 
as a Schechter function 

(j)(L2so)dL = <j>*(L/L*) a exp-L/L*dL . (1) 

We chose the restframe UV luminosity at 280 nm, L280, since it 
can be used as a proxy for the u nobscured star-formation rate using 
the relation o f lKennicuttl fl998) 

SFRtMoyr" 1 ) = 1.4 x 10~ 28 L 28 o/(erg s" 1 Hz" 1 ) . (2) 

Evidently, the total star-formation rate of a galaxy depends also on 
the amount of obscured star formation, which could be measured 
from the thermal FIR dust emission or crudely estimated using UV 
spectral slopes. 

However, both CC SNe and LGRB afterglows are subject to 
host galaxy extinction as well, just as expected for events occur- 
ing in very young stellar populations. In fact, star-bursting galaxies 
do not only have an above average proportion of obscured star for- 
mation, but also an increased rate of obscured, o ptically invisible 
CC S Ne, which can be found by NIR monitoring dMannucci et al.l 
2003). Such SNe are absent in the sample used here, which is 
based on a restframe UV-optical search. Present data for LGRB 
afterglows show no evidence for high levels of dust obscuration, 
while a fraction of optically dark LGRBs may still be due to high 
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Table 1. UV luminosity functions of galaxies at z : 
and GEMS (transformed from h to hn)- 



: 0.7 from COMBO-17 



sample 


M* 80 - 5 k 




a 


all galaxies 


-20.45 ± 


0.16 


-0.75 ±0.20 


irregular galaxies 


-20.11 ± 


0.17 


-1.31 ±0.20 



dust extinction in the host galaxy dKlose et al.ll2003l) . However, the 
host galaxies of known LGRBs show no signi ficant signs of dust- 
obscured star formation dLe Floc'h et al J 20061) . 

Hence, we conclude that a comparison of UV-based star for- 
mation measures with UV-optical detection rates of stellar explo- 
sions is fair. It would be next to impossible to attempt an accu- 
rate extinction correction, which would need to take into account 
changes in the average extinction with progenitor lifetime as well 
as detailed detection efficiencies for events of different extinction 
levels. 

For our model, we first calculate the distribution of UV lumi- 
nosity density j280 cumulative over UV host luminosity L2&0 using 
the integral 



J280(i2 



1/280 0(^280) d-I/280 



(3) 



Given the Kennicutt law, this quantity is an estimate of the unob- 
scured star-formation density in hosts of L280 > £<280,iim- Hence, 
it also describes the host distribution of any explosive event which 
is proportional to the star-formation rate. As an example, CC SNe 
are broadly expected to be an unbiased tracer of star formation, and 
should have a host distribution described by ^3). 

In a second step, we transform the UV luminosity axis L28O 
into Lv to match the given luminosity data of the host sample. 
For this purpose, we need to know the (M280 — My )-colours 
of star-forming galaxies. These data have been obtained by the 
COMBO-17 surve y and are publicl y available for the Chandra 
Deep Field South jWolf et alj|2004l) covering the entire GEMS 
field. The galaxy sample used for the above LFs can be split into 
a red sequence of mostly passively evolving galaxies and a blue 
cloud of mostly star-forming galaxies. For star-forming galaxies at 
2 ~ 0.7, the colour-magnitude relation is (see Fig. 5 of W05) 



Ah 



My = -0.11(M 28 o - 5 log h ro + 20) 



(4) 



Hence, the lower integration limit in the calculation of the luminos- 
ity density can be substituted as 



M v 



= Mo 



+ 0.11(M 2; 



5 log h 70 + 20) 



(5) 



This is equivalent to stretching the scale of the luminosity axis by 
11% and holding it fixed at My = M280 = —20, where the aver- 
age colour of star-forming galaxies is zero. 



3.2 Predictions by galaxy type 

It has often been reported, most recently by F06, that CC SNe ap- 
pear in all types of star-forming galaxies, while the afterglows of 
LGRBs are almost exclusively found in morphologically irregular, 
faint blue galaxies. Although galaxy morphology may not be a suit- 
able physical factor for the evolution of high-mass stars in a galaxy, 
it may be correlated with physical factors: irregulars are mostly of 
lower-luminosity and hence tend to be of lower metallicity, which 
may well be the physical cause of the increased LGRB rate. 
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Figure 1. Cumulative luminosity distributions of CC SN and LGRB 
host galaxies (step functions), together with distributions for the UV- 
derived star-formation density in all galaxies and in irregular galaxies alone 
(curves). 



Table 2. L-Z Q -relations at z ~ 0.7 from Kobulnicky & Kewley (2004) 
using a functional form of 12 + log (O /H) = a + b Mb ■ 



fit 


a 


slope b 


A 


6.37 


-0.117 


B (best) 


3.88 


-0.239 


C 


1.37 


-0.363 



We now briefly try to confirm the picture suggested by F06, 
and test a toy model in which CC SNe arise unbiased from (un- 
obscured) star formation in all galaxies while LGRBs arise unbi- 
ased from irregular galaxies. We model the host luminosity distri- 
butions using the parameters of the LFs obtained by W05 for ir- 
regular galaxies and all galaxies (see Tab. Q}. Note that the model 
is only a normalised cumulative luminosity distribution. Hence, it 
does not constrain overall rates, and it also does not depend on the 
normalisation <f>* of the LF, but only on the slope a of the faint end 
and the location M* of the break in the LF. 

This toy model is compared with the sample data in Fig.Q] We 
find that CC SNe match the prediction extremely well and conclude 
that they can be considered an unbiased proxy for star formation in 
all galaxies within the limits of this small data set. Also, the pre- 
diction for irregular galaxies matches the host data for LGRBs rea- 
sonably well, although there is some deviation at low luminosities. 
Formally, both comparisons are entirely acceptable: a Kolmogorov- 
Smirnov test yields values of Dks = 0.12 and 0.16, respectively 
(corresponding to rejection probabilities of only 4% and 22%). 

We note, that the LFs are well-determined at My < —18.5 
but extrapolated with a Schechter function at fainter magnitudes. 
Any errors in this extrapolation will translate into a simple rescal- 
ing of the cumulative host prediction at My brighter than —18.5, 
whereas the fainter section would change its shape. We choose to 
ignore the errors in the previous qualitative analysis and present 
a more quantitative study, using metallicity as the physical factor 
controlling the LGRB efficiency, in the following sections. 
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3.3 Luminosity-metallicity relations and dispersion in the 
stellar population 

In this section, we refine our toy model by considering as progen- 
itors only sub-populations selected by metallicity. Our use of the 
term metallicity denotes the abundance of oxygen in the ionised gas 
phase as a tracer of metallicity in young high-mass stars. We use 
the O/H ratio measured via the i?23-index on the KK04 calibration. 
This index is determined from flux ratios in the prominent ionized 
oxygen and H/3 emission lines of star-forming galaxies, which are 
relatively easy to observe even in faint galaxies at higher redshifts. 
Of course, it would also be attractive to look at abundances of other 
elements. However, these are more challenging to measure, and we 
believe that, with the present small sample of objects, we would not 
be able to improve our study significantly. As an abbreviation we 
will use the notation 



Zo = 12 + log (O/H) 



(6) 



Recent determinations of the luminosity-metallicity (L-Zo) or the 
stellar mass-metallicity (M« -Zq) relation include the work by 
iKobulnickv & Kewlevl ( 12004 herafter KK 041 on the basis of the 
GOODS survey, as well as lSavaglio et all d2005h using the GDDS 
survey. We opt for using the L-Zo -relation, which relates more di- 
rectly to our data and model. For star-forming galaxies at z ~ 0.7, 
the COMBO-17 data suggest on average Mb — My, and we 
use the KK04 relations without further modifications. Their rela- 
tions are determined in redshift slices, and conveniently their slice 
z = [0.6, 0.8] corresponds precisely to the redshift range in which 
the LF was determined by W05. Altogether, three linear fits are 
shown in their Fig. 11: a Zo = /(A/ s )-fit (A), an Mb = / (Zo)- 
fit (C) and a bisector fit (B), which can all be expressed as 



Z Rt = 12 + log (O/H) = a + b x Mb 



(7) 



The bisector fit has a slope of b = 0.239, between the two one- 
sided fits, and is presumably the most realistic one of the three. We 
also use the two one-sided fits to define a confidence interval for 
the relation (see Tab.[2]for fit parameters). 

A recent determination of the L-Zo -relation in local dwarf 
galaxies suggests a slope of b ~ 0.30 e xtending across the lumi- 
nosity range from M B = -11 to -19 jLee et alj|200^) . At fixed 
luminosity, the metallicity has a la-scatter of 0. 16 dex, the same as 
that determined in a l arge sample of luminous galaxies from SDSS 
dTremonti et ali20 04». The scatter of the mass-metallicity relation 
is even less (0.10), and in the absence of any higher-redshift data 
we assume that it is well-behaved also for dwarfs at z ~ 0.7. In 
any case, both the L-Zo -relation and the UV luminosity function 
are just extrapolated below My — —18.5 in our analysis. 

Below we first investigate the effects of a sharp metallicity 
cutoff for the production of LGRBs in our model. If we ignored 
both the scatter in the L-Zo -relation and the internal variations of 
Zo among the young stellar population within any given galaxy, a 
sharp Zo -cutoff would cause a sharp My -cutoff in the host galaxy 
sample. While we do not have perfect knowledge of the overall Zo 
scatter, it certainly needs to be taken into account for the prediction 
of cumulative host luminosity distributions in the presence of in- 
dependent event efficiencies. 

Even for the same data, the scatter of the L-Zo -relation 
for average galaxy metallicities varies with the slope of the fit. 
Roughly, it is around ±0.3 dex end-to-end for a given luminosity at 
z ~ 0.7, consistent with a lcr-scatter of 0. 16 at low redshift. As far 
as the Zo -spread of star-forming gas within a galaxy is concerned, 



we lack useful knowledge at z ~ 0.7 owing to the limited spatial 
resolution of spectroscopic data, especially for smaller galaxies. 

Instead, we look at the present-day Magellanic clouds for 
clues on the Zo variation within irregulars. The LMC appears 
to resemble a typical LGRB host at z ~ 0.7 given its lumi- 
nosity of My ~ —18.5. The LMC and SMC have measured 
metallicities in young stars and the gas phase of 0.5 solar and 
0.2 solar, r espectively, and both h ave internal la-variations of 
< 0.2 dex ( lRussell&Besselllll989l) . Dispersions in Milky Way- 
size disk galaxies are not expected to be larger either, given efficient 
gas mixing. 

In principle, we need to convolve the distribution of the inter- 
nal Zo variation and the scatter in the L-Zo -relation to obtain a 
realistic cosmic p(Za\Mv) -distribution. But in the absence of any 
accurate values for the former, we consider two alternative scenar- 
ios for the combined metallicity variation by using the form 



dz{Zo,M v ) 
p{Z \Mv) 



Zo - Z flt (M v 



if|dz|>l, 
cosdz if \dz\ ^ 1. 



(8) 
(9) 



where Zat is taken from the L-Zo -relation and ±w is the end- 
to-end spread in the variation applied at any fixed My. We opted to 
use a cosine function, because it is trivial to integrate analytically. 
Two other simple alternatives appear unrealistic, a box function be- 
cause of its steep edges and a Gaussian because of its infinitely 
far-reaching wings. 

Overall, we consider w = 0.4 dex a realistic value for the 
combined spread at fixed host luminosity, which corresponds to 
about ±0.25 dex for the FWHM of the distribution p. Just for illus- 
tration we also explore an extreme value of w = 1, which translates 
into an extreme and unrealistic metallicity spread from ten-fold be- 
low to ten-fold above the fit value at any given luminosity. This will 
allow even very luminous galaxies in our model to form a fraction 
of very low-Zo stars. 



3.4 Event efficiencies and their dependence on metallicity 

Here, we specify our descriptions for the dependence of the LGRB 
efficiency on the metallicity of young stars. We consider two simple 
one-parameter models and one two-parameter model, which is a 
compromise between the first two. Any adhoc model that is more 
complicated than these makes little sense given the small number 
of 16 LGRBs in the sample. Our adopted prescriptions are: 

• a sharp high-Zo cutoff: 



• a power law, which unfortunately diverges at low Zo : 
lo g7? (Z ) = -cZo. (ID 

• a broken power law, which is flat as in J 10b below a break 
metallicity Z* and declines as in Jilt above the break: 



logTj(Zo) = 







if Zo - Z* < 0, 



-c(Z -Z„) ifZ o -Z»>0. 



(12) 



Note, that for the single power law no normalisation is required 
as we are not dealing with absolute LGRB rates but only host dis- 
tributions. Also, model d 12b contains the other models as special 



© 0000 RAS, MNRAS 000, 000-000 



The metallicity dependence of the LGRB rate 5 




-22 -20 -18 -16 -22 -20 -18 -16 

M v -5 log h 70 M v -5 log h 70 



Figure 2. Top left: Cutoff models for a sequence of metallicity cutoffs at (Z\ iBl = {9.0, 8.7, 8.4, 8.1}). The reference model for all galaxies has no metallicity 
dependence. Bottom left: Power-law efficiencies with c = {—0.2, —0.4, —0.6, —0.8}. All models predict too many bursts in faint galaxies. This is a 
consequence of the divergence to infinity of the LGRB efficiency at low metallicity. Varying w has no effect on power-law models. Right: Broken power 
laws: The softer the break in the efficiency function, the lower the break metallicity in the best-fitting models. A hard cutoff (c = —3, top) is best fit with 
Z\x m ~ 8.7 (i.e. solar), while a very soft break (c = — 1, bottom) is best matched at Zi; m ~ 8.2. In the latter case the LGRB efficiency at solar metallicity is 
still a third of the plateau value. 



cases, because it converges into JlOt for c — > — oo, and into dl It 
for — > — oo. 

The prescriptions for event efficiencies rj{Zo) and for the 
metallicity distributions of star-forming gas p(Zo\Mv) in hosts of 
given luminosity My are then included into our prediction of the 
V-band host luminosity distribution, using the new integral 

j28o(Lv,lim) = dI/28o/d£vX 

/ / p(Zo\L v ) v (Zo)L 2 s (L v )4>(Lv)dL v dZo ■ 

"-Lvlim — 00 

Again, this figure is now proportional to the LGRB event rate ex- 
pected for the assumed efficiency function. 



4 COMPARISON WITH THE DATA 

In this section, we compare our metallicity-dependent models with 
the host galaxy data. However, in Sect. 5.2 we will discuss empir- 
ical evidence that the LGRB hosts discussed here are on the L-Z 
relation as far as we know, supporting the validity of our basic ap- 
proach. 



4.1 Cutoff models 

Here, we use a constant efficiency up to a cutoff metallicity defined 
in eq. UOi and consider a sequence of cutoff values of Z\i m = 
{9.0, 8.7, 8.4, 8.1}, i.e. roughly twice solar, solar, half solar and a 
quarter solar. The top left panel in Fig.|2]shows the resulting models 
alongside the reference model with all galaxies that fits the CC SN 
hosts. The most likely cutoff metallicity can be readily obtained 
from the plot as ~ 8.7. Already a cutoff at 8.4 clearly underpredicts 
the host luminosities. 

We also investigated the dependence of this results on uncer- 
tainties in the L-Zo -relation by using the limiting fits A and C 
from Tab. [2] to define fiducial confidence intervals. We find that 
changes in the slope of the L-Zq -relation translate directly into 
changes in the luminosity spread of the curves predicted for a fixed 
set of cutoff metallicities, which is expected. Fortuitously, the data 
lie in a range that appears to be not very sensitive to the choice of 
the L-Zo-fit. This is a result of the pivotal point in the L-Zo-fit 
being close to M = —20. There are not enough GRB hosts far 
away from the pivotal point of the L-Zo-fit to let the slope make 
a large difference for the best-fitting cutoff metallicity. Instead, the 
slope of the relation affects the formal errors on the cutoff, which 
one obtains from a fit to the data (see Sect. 4.4). 
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4.2 Power-law efficiencies 

As an alternative to a sharp cutoff model, we now test an efficiency 
function that smoothly changes with metallicity and choose a sim- 
ple power law. The reference model of a constant efficiency for all 
galaxies is always shown for comparison, and corresponds to c = 
in eq. i ll It . Fig.|2]further shows the curves for power-law slopes of 
—0.2, —0.4, —0.6 and —0.8, which deviate increasingly from the 
reference model as high-Zo GRBs are more and more suppressed 
relative to low-Zo events. 

Any variation in w has no effect on the curves, because the 
symmetric form of the Zo distribution is averaged out with a 
power-law efficiency. Hence, there is no need to reproduce the 
curves for any other w. Furthermore, the slopes c of the efficiency 
curves are degenerate with the slope of the L-Zo -relation. So, fits 
A or C will produce the same curves, if suitable c-values are cho- 
sen. All curves for power laws predict a substantial contribution 
from extremely low-luminosity galaxies, which are not observed 
among GRB hosts. This is due to the unrealistic divergence of the 
GRB efficiency towards zero metallicity. 



4.3 Broken power laws 

The broken power law with flat \ow-Zq efficiency combines the 
properties of the two previous models, a gradual decline towards 
higher metallicity and a constant, non-diverging efficiency at lower 
metallicity. Hence, this model avoids the unrealistic and problem- 
atic divergence of the single power law, while being less restrictive 
than the cutoff model, providing the "best of both worlds". 

The model is now described by two parameters, which pro- 
vides more freedom to fit the data, but also introduces some degen- 
eracy. While hard breaks are very similar to the cutoff models, the 
metallicity break locations for softer breaks become increasingly 
less constrained. Especially, with only 16 GRBs the data do not 
provide sufficient constraints for the model. Conversely, we con- 
clude that these data do not yet rule out a variety of models. 

Fig. [2] shows results for two broken power laws to illustrate 
the trends with the slope c. The sequence of models shows breaks 
at = {9.0, 8.7, 8.4, 8.1}; a constant efficiency matching the CC 
SN data is also shown. We find that the data can not differentiate 
between a soft and a hard break. Instead, a hardening of the break 
only leads to higher break metallicities. While the assumption of 
a hard cutoff suggests Z* ~ 8.7, softening the break to c = —1 
moves Z* to ~ 8.2. 



4.4 Kolmogorov-Smirnov tests and results at redshift < 1 

In the present analysis systematic uncertainties appear to dominate, 
ranging from the L-Zq -relation over the typical metallicity spread 
in the star-forming gas of galaxies at 2 ~ 0.7 to the form of the 
faint end of the ultraviolet galaxy luminosity function. However, 
it is straightforward to obtain quantitative confidence intervals on 
Zjin, (as shown in Fig. [3} or Z* (Fig. [4]l for a particular model, 
using Kolmogorov-Smirnov (KS) tests. We repeat, that the KS test 
yielded a highly acceptable value of Dks = 0.12 for the match 
between the CC SNe and the reference model of 'all' galaxies. 

Fig.[3]shows the results for the cutoff models. The solid curves 
show the three different L-Zo -relations with a realistic metallicity 
spread of w — 0.4. All three favour a cutoff near solar metallicity 
and differ mostly in the width of the confidence interval, which is 
a result of the different slopes of the relations. We obtain Zu m = 



id 




Figure 3. Kolmogorov-Smirnov tests: Cutoff models with w = 0.4 (solid 
curves) and different L-Zo-fits. The dashed curve depicts a model with 
extreme metallicity inhomogeneities w = 1. 
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8.7±0.22 at 95% c.l. for the most likely L-Z -relation B or ±0.3 
if allowing for any relation. 

But how would these results be affected if the internal metal- 
licity dispersion of star-forming galaxies was truely much larger 
than ever anticipated? Would this allow the low-metallicity tails 
in luminous galaxies to produce GRBs and make the host data 
consistent with much lower metallicity cutoffs? For illustration, 
we modelled an extremely wide metallicity spread of the star- 
forming populations at given host My, from l/10th of the mean 
to lOx the mean. While such a model is unrealistic, the loca- 
tion of the cutoff metalllicity is only lowered by ~ 0.25 dex (see 
dashed line in Fig. [3}. Invoking such extreme metallicity inhomo- 
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geneities among young stellar populations still rules out cutoffs of 
Zn m < 8.1 « Zq/4 at 95% confidence. 

Fig-Sillustrates the broken power-law models. Softening the 
break, from a hard cutoff (c = — oo) over c = —3 to a slope of 
c = — 1 moves the favoured break metallicity to lower values and 
widens the confidence interval. A more gradual change in GRB 
efficiency leads to a more gradual rise of the cumulative host lumi- 
nosity distribution. A soft break is also a less pronounced feature in 
the efficiency function rj(Zo). The best-fitting 77 functions shown 
in the bottom panel illustrate that even softer breaks require the 
r,(Z Q ) « 1/3. 

In summary, we have explored for the z = [0.2, 1.0] sample 
the impact of the two most important uncertainties by comparing 
results for different sets of assumptions. Three alternative L-Zq- 
relations and different degrees of internal metallicity scatter in the 
star-forming gas are put in contrast in Fig. [3] Varying the L-Zo- 
relations around the pivotal points provided by the observations of 
z ~ 0.7 galaxies appears to have very little effect. Changing the 
internal metallicity scatter of typical galaxies can have a significant 
effect when assuming unconventionally extreme degrees. 

We further calculated the fractions of the young stellar popu- 
lation that can produce LGRB progenitors for z ~ 0.7, the mean 
redshift of our LGRB sample, given the various efficiency pre- 
scriptions. In the metallicity-independent model where all galax- 
ies contribute and match the CC SN distribution, this fraction is 
defined as equal to 1. In the cutoff model the best-fitting value 
Z\i m = 8.7 ± 0.3 corresponds to fractions of 0.55 ± 0.25. This 
means that because of the metallicity constraint the LGRB produc- 
tion efficiency is reduced by a factor of about one half compared to 
the efficiency if there were no metallicity constraint. This number 
is also consistent with half of the star formation and half of the CC 
SNe originating from galaxies with lower metallicity than solar, or 
with lower luminosity than Mb = —20. In the broken power-law 
models we find very similar fractions independent of the ambiguity 
between break softness and break location. In the single power-law 
model the fractions can not be normalised because the efficiency 
diverges at low Zq . However, these models did not fit the data for 
the same reason. 
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Figure 5. Top panel: Comparison of the fiducial isotropic LGRB energy 
i?i so with the host metallicity, as estimated from the L-Zq -relation for 
z > 0.2-bursts (black points) and as measured from spectroscopy for four 
z < 0.2-bursts (grey points). Bottom pa nel: Comparison of the LGRB 
energy E corr corrected for beaming by Ghirlanda, G hisellini & Lazzatil 
120041 with the host metallicity. 



4.5 High-redshift host galaxies 

The purpose of restricting the GRB sample to low redshift was 
to compare it directly against the CC SN sample and to quanti- 
tatively investigate the role of metallicity in a regime where we 
have ample knowledge about the star formation and metallicity 
of the galaxy population. After obtaining a possible value for 
a metallicity cutoff around Zn m ~ 8.7, we can check qualita- 
tively whether this is in agreement with the observations of higher- 
redshift hosts. The F06 sample contains eight host galaxies in the 
range z — [1.5, 3.5], half of which have luminosities brighter than 
Mv = —20. This reflects the fact, that at higher redshift larger 
galaxies wer e more actively forming s tars than they do at lower- 

redshift (e.g. lPerez-Gonzalez et alj2005|). 

A recent study of z ~ 2-galaxies terb et alj2 006) showed that 
there is a mass-metallicity relation in place but no well-defined L- 
Zo -relation due to strong variations in the M/L-ratios of galaxies. 
In their study, the most massive galaxies reach Zo = 8.6, which 
is on a T e -based calibration and probably equivalent to ~ 8.9 on 
the KK04 scale. A limit at Zo ^8.7 would correspond to a mass 
limit of M»/M < 10 10 ' 5 . In fact, the sample contains plenty of 
luminous galaxies below this mass with My up to —22. Hence, 
our findings from the z < 1-sample are not in conflict with the 



luminosities of higher-redshift hosts. On the other hand, the Erb et 
al. sample contains galaxies with metallicities above our estimated 
cutoff, suggesting that at z ~ 2 LGRBs are not entirely unbiased 
tracers of star formation, which they might be at yet higher redshift. 



5 DISCUSSION 
5.1 Caveats 

One difficulty with the present analysis could be host confusion. If 
the GRBs took preferentially place in low-luminosity companions 
to the formally identified hosts, then the observed host luminos- 
ity distribution would be biased to higher luminosities, suggesting 
higher metallicity cutoffs. It is beyond the scope of this paper to 
investigate whether HST observations could resolve such dwarfs 
from their larger, mostly LMC-type companions at z < 1. How- 
ever, if most LGRBs occurred in such small galaxies, one would 
also expect a large population of sufficiently isolated small dwarfs, 
where no confusion would take place. But the small number of ob- 
served isolated hosts with Mv > — 17 leaves little room for any 
adjustment of our results. 

On the contrary, it is possible that more massive LGRB hosts 
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are missing from the sample: luminous dust-obscured star-bursting 
galaxies harbour a significant fraction of the overall cosmic star 
formation. If LGRBs occur in these, their optical afterglows may 
be invisible, which would depopulate the host sample specifically 
at the bright end. The CC SNe host sample would be less affected 
by this effect, as CC SNe have a lower average progenitor mass 
and hence longer average time delay from formation to explo- 
sion. This allows the progenitors to leave their immediate forma- 
tion environment and travel into less obscured regions. The effect 
of progenitor mass on the distance to the likely formation site has 
been shown to apply in t he comparison of SNe II and SNe Ib/c 
jjames &~A nderson 2006). The sample analysed here lacks by def- 
inition all so-called 'dark bursts', which may be related to massive 
dusty galaxies: GRB 020127 only produced an X-ray afterglow, 
and the host was identified as a z ~ 1.9 dusty ERO galaxy with 
5 L* luminosity teergeretalj|2006h .l Levan et al. (2006) observed 
an extremely red (R — K ~ 6) afterglow to GRB 030115 in an 
ERO host galaxy with R — K ~ 5, suggesting that this burst was 
highly obscured by dust. The afterglow was only identified with 
NIR observations and would have been missed by optical follow- 
up alone. Again, for the more recent GRB 050223 only an X-ray 
afterglow was seen by Swift, but no optical afterg low. The host wa s 
shown to be a dusty galaxy with Av > 2 jPellizza et al.l l2006). 
If there is a significant observational bias against LGRBs in mas- 
sive dust-enshrouded galaxies, then the metallicity-dependence of 
LGRB rates may be rather weaker than our findings presently sug- 
gest. 

Uncertainties also exist with respect to the calibration of 
oxygen gas-phase metallicities and on how the oxygen abun- 
dance is representative of other elements. E.g., a-elements are 
often found to be enhanced in systems where the star-formation 
timescale is much shorter than the timescale for Type la super- 
novae. Depending on the progenitor model for LGRBs, different 
elements might be most responsible for affecting opacities in stel- 
lar atmospheres/winds and hence stellar evolution. However, a- 
enhancements are typically observed in old ellipticals which are 
not at all common LGRB hosts or typical star-forming galaxies, 
and thus should not affect our analysis appreciably. 

In summary, solar metallicity appears to mark a pivotal point 
of roughly constant efficiency apparently independent of the soft- 
ness of the break. At z ~ 0.7, this metallicity seems typical for 
galaxies of Mb — 5 log hro ~ —20 in either L-Zo -relation. The 
number ratio of CC SNe to either side of this host luminosity cut is 
1:1, whereas it is 1:3 for LGRBs. If there was any doubt about the 
calibration of the oxygen abundance at 2 ~ 0.7, or the relevance 
of oxygen in comparison to other elements, then we could rephrase 
our main conclusion such that a most likely cutoff for LGRBs is 
around the mean metallicity (by any measurement) of Mb — — 20 
galaxies at z ~ 0.7. 

5.2 Are GRB hosts on the luminosity-metallicity relation? 

Several recent studies have aimed at constraining LGRB 
progenitors from dir e ct measurements of host metallicities 
(Solle rmanetai]|2005l : ISavaglio, Glazebrook & Le Borgnel 120061: 
IStanek et alj |2006|)~ using very small samples: ISollerman et al.1 
d2005h present measur ed metallicities for three LGRB hos t s with 
values of ~ 8.2 . . . 8.7. ISavaglio. Glazebrook & Le Borgnel (2006) 
consider seven hosts at z — [0.4, 1.0] and find mean galaxy metal- 
licities of 8.3 to 8.55 for five of them, while two hosts are clearly 
super-solar. They furthermore claim, that the hosts fall right onto 
the regular mass-metallicity relation of normal star-forming galax- 



ies. This host sample with measured metallicities lacks the lowest- 
luminosity hosts and so contains preferentially higher metallicity 
galaxies. Hence, their median metallicity is not representative, but 
the statement on consistency with the mass-metallicity relation is 
highly relevan t. 

ISavaglioi (2006) give a number of host metallicity determina- 
tions, which are collected from the literature and converted onto the 
same (KK04) calibrator, which is also the calibrator used in this pa- 
per for the L-Zo -relation. Five of the above hosts are actually part 
of the LGRB host sample used here. Their spectroscopic metallici- 
ties are compared to the estimated values we have assigned to them 
using L-Zo-fit B. We find them all to be within ±1/3 dex of the 
relation with a mean of —0.034 dex, which is consistent with no 
bias within the statistical power of five objects. 

This results is not surprising. A physical model for GRB pro- 
genitors may incorporate an explicit metallicity dependence, but 
we can not expect the progenitor to know explicitely about galaxy 
parameters such as luminosity or morphology. At a fixed metallic- 
ity level, we would then expect bursts predominantly from those 
galaxy mass bins containing the highest star-formation rate. Here, 
both the shape of the mass function and the SFR trends with mass 
play a role. Presumably, the declining mass function and the in- 
creasing SFR with mass cancel to some extent. KK04 searched 
for trends of galaxy properties within the scatter of their L-Zq- 
relation. While they found no such trends, they particularly did 
NOT find strongly star-forming galaxies on the bright or metal-poor 
side of the relation. However, among very low-mass dwarf galax- 
ies it is conceivable, that strong variations in SFR with time lead to 
such strong variations in mass-to-light ratio, such that CC SNe and 
LGRBs are found predominantly on the bright side of the relation. 
Presumably, this would have no effect for our analysis, as the mean 
metallicity bias would be small and remain within the flat portion 
of the parametric efficiency functions considered in this paper. 



5.3 Global vs. local metallicity measurements 

We would like to point out, that even spectroscopic, but spatially 
unresolved, measurements of a mean host metallicity do not di- 
rectly reflect the metallicity of the young stellar population near 
the LGRB progenitor. In fact, the error in the progenitor metallicity 
is dominated by the internal Zo dispersion of the galaxy, which is 
estimated to be as large as the dispersion of galaxy-averaged Zo 
values in the L-Zo-relation (< 0.2 dex). 

If we make no prior assumptions on LGRB hosts and accept 
the result that LGRB hosts follow the usual L-Zo -relation, then 
we can in turn estimate even inidividual host metallicities from this 
relation. These individual estimates then have errors on the order 
of the dispersion in the L-Zo -relation, i-e. < 0.2 dex. 

In other words, owing to the internal metallicity dispersion, 
our indirect estimate of the progenitor metallicity via the L-Zq- 
relation of the host galaxy should be almost as useful as an un- 
resolved spectroscopic observation, with an error that is roughly 
larger by a factor of V2. If the internal metallicity dispersions ex- 
ceeded 0.2 dex, then the determination of integrated metallicities 
of individual LGRB hosts would have very little value, because it 
could as well be estimated from an L-Zo -relation, while the error 
was mostly internal in origin. 

Therefore, major progress can mostly be expected from spa- 
tially resolved metallicity measurements that focus on the imme- 
diate burst environment. However, such observations are currently 
not possible for GRBs at cosmological distances. 
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5.4 A metallicity-LGRB energy relation? 

IStaneketalJd2006h suggest a tentative relation between host metal- 
licity Zo and the isotropic energy associated with the LGRB Ei so , 
which they interpret as a metallicity cutoff around 0.15 solar for 
regular cosmological bursts. This conjecture is based on one ob- 
ject out of five, where a high LGRB energy coincides with low 
host metallicity. While this object was assigned a much higher host 
metallicity in a previous analysis (: Solle rman et al. 2005). the mass 
of the galaxy makes that interpretation very unlikely. While the 
claimed relation relies on the significance of this single object, it 
is unclear whether the other (underluminous) LGRBs are represen- 
tative of the more energetic and distant cosmological bursts. GRB 
060218, e.g., could not even have been detected by Swift at redshifts 
2 > 0.05. 

Even in the absence of spectroscopic metallicity measure- 
ments for the individual hosts, we can tentatively investigate such a 
relation for the full sample with known i5i so and host luminosities 
using the L-Zo -relation. Fig.[5]shows the result for 13 GRBs from 
our sa mple (black data points) which have Ei so values in the liter- 
aturejAmati 2006). We also add the five local objects discussed by 
Stane k et al .(2006) based on their actual host Zo-data as grey data 
points. At least for the cosmological bursts (black) no relation is 
apparent. However, due to the expected significant variations in jet 
geometry and viewing angle among the bursts, we would not expect 
the fiducial isotrop ic energy estimate to show very c l ear tre nds with 
other parameters. Ghirlanda, Ghiselli ni & Lazzatil d2004h present 
LGRB energies -E C orr that are corrected for jet beaming effects. 
Correspondingly, their Ghirlanda relation E corr vs. Speak is much 
more clearly defin ed than the original Amati relation i5i so vs. Speak 
dAmati et alj2002t) . where -E pC ak is the restframe photon energy at 
the peak in the vF v spectrum. In the bottom panel of Fig. [5] we 
show -E cc> rr vs. the host metallicity for the bursts with available 
data. We find no trend in these data. 



5.5 Theoretical implications 

From a theoretical point of view, it is not surprising that metal- 
licity plays an important role for LGRB progenitors. In one of the 
most promising progenitor scenarios, the collapsar model (Woosley 
1993; MacFadyen & Woosley 1999), the progenitor is the rapidly 
rotating core of a massive star. On the other hand, the core of a mas- 
sive star is believed to efficiently lose angular momentum during 
its evolution: both within the star by hydrodynamical (e.g. Heger, 
Langer & Woosley 2000) and magnetohydrodynamical processes 
(e.g. Spruit 2002) and subsequently from the surface of the star in 
the form of a stellar wind. As a consequence, it seems unlikely that 
most massive stars will still have cores at the time of core collapse 
that rotate as rapidly as is required in the collapsar model (Heger, 
Woosley & Spruit 2005). 

Indeed, this is also not necessary, since it is clear that LGRBs 
are rare events, associated with less than 1 in 1000 core-collapse su- 
pernovae (Podsiadlowski et al. 2004b). One of the key unresolved 
questions is what are the special circumstances that produce LGRB 
progenitors: are these special conditions in single stars (e.g. Yoon & 
Langer 2005; Woosley & Heger 2006) or do they require particular 
binary channels (e.g. Fryer & Woosley 1998; Izzard, Ramirez-Ruiz 
& Tout 2003; Fryer & Heger 2005; Petrovic et al. 2005; Detmers 
et al. 2006; Fryer et al. 2006; Podsiadlowski et al. 2006)? One of 
the major problems in many of these models is that mass loss, both 
in the red-supergiant phase and, in particular, in the Wolf-Rayet 
phase, leads to very efficient angular-momentum loss. Since these 



wind mass-loss rates are dependent on the metallicity (typically 
M oc Z°- 5 ~°- 7 \ e.g. Vink & de Koter 2005), massive stars with 
lower metallicity are more likely to have rapidly rotating cores at 
the time of core collapse. This is a generic advantage for many of 
the proposed models, both single and binary. 

Specifically, Yoon & Langer (2005) and Woosley & Heger 
(2006) recently proposed that a low-metallicity, very rapidly rotat- 
ing star may evolve essentially homogeneously during its early evo- 
lutionary phases and may avoid a red-supergiant phase altogether, 
in which most of the angular-momentum loss from the core tends 
to take place. As a consequence, these stars would still preserve 
rapidly rotating cores at the time of collapse, fulfilling one of the 
key conditions in the collapsar model. Similarly, in some of the 
most promising binary scenarios, in which two massive stars (al- 
most) merge or interact tidally, it is advantageous or even necessary 
that this occurs late in the evolution of one of the stars, i.e. after he- 
lium core burning (so-called case C mass transfer). The range of 
masses for which case C mass transfer occurs and allows the for- 
mation of a black hole is a strong, non-linear function of metallicity 
(Justham et al. 2006), again favouring low-metallicity progenitors. 

As this discussion illustrates, the dependence of the 
LGRB rate on host metallicity can provide an important con- 
straint on the various p r opos ed proge nitor models. Indeed, 
Hirschi, Meynet & Maede j (2005) and I Yoon. Langer & Normanl 
d2006h have estimated that, when the effects of magnetic fields 
are included, the single, low-metallicity, hi gh-rotation model 
(Yoon & Langer 2005; Woosley & Heger 2006) requires a metal- 
licity less than l/5th solar. This already appears inconsistent with 
the constraints derived in this paper, which seem to rule out any 
models that require a metallicity significantly less than 1/2 solar. 

The ongoing metallicity calibration debate might render our 
limit to lie ~ 0.3 dex lower, if T e -based methods were correct. 
There are now explorations of a third method for measuring metal- 
licities, the O Hrl method dPeimbert etall kooe), which is aimed 
at overcoming the shortcomings of both previous methods. Early 
results indicate that this method gives values in between the two 
debated versions, but closer to the R23 method. This would leave 
our results largely unchanged. 
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